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q ■ Abstract: 

I discuss possibilities to observe the instanton-induced contributions to deep inelastic 
scattering which correspond to nonperturbative exponential corrections to the coefficient 
functions in front of parton distributions of the leading twist. 
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1. The deep inelastic lepton-hadron scattering at large momentum transfers Q 2 
and not too small values of the Bjorken scaling variable x = Q 2 /2pq is studied in much 
detail and presents a classical example for the application of perturbative QCD. The 
celebrated factorization theorems allow one to separate the Q 2 dependence of the structure 
functions in coefficient functions Ci(x, Q 2 / /i 2 , a s ([i 2 )) in front of parton (quark and gluon) 
distributions of leading twist Pi(x, /i 2 , a s (fi 2 )) 

F 2 {x, Q 2 ) = EiQix, Q 2 /fx 2 , a s (/i 2 )) <g> P(x, /i 2 , a s (fx 2 )), (1) 

where 

C{x)9P(x)= t ^-C(x/y)P(y), (2) 
Jx y 

the summation goes over all species of partons, and \i is the scale, separating "hard" and 
"soft" contributions to the cross section. At fi 2 = Q 2 the coefficient functions can be 
calculated perturbatively and are expanded in power series in the strong coupling 

C(x, 1, a s (Q 2 )) = C (x) + ^^-C x {x) + (^1) ' C 2 (x) + ... (3) 

7T y 7T J 

whereas their evolution with fj, 2 is given by famous Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi equations. Going over to a low normalization point fi 2 ~ lGeV, one obtains the 
structure functions expressed in terms of the parton distributions in the nucleon at this 
reference scale. The parton distributions absorb all the information about the dynamics 
of large distances and are fundamental quantities extracted from the experiment. Pro- 
vided the parton distributions are known, all the dependence of the structure functions 
on the momentum transfer is calculable and is contained in the coefficient functions Cj. 
Corrections to this simple picture come within perturbation theory from the parton distri- 
butions of higher twists and are suppressed by powers of the large momentum Q 2 . These 
higher-twist contributions are relatively well understood theoretically, and, unfortunately, 
very poorly experimentally. I have nothing to say about them in this talk. 

The picture described above presents a part of the common wisdom about hard pro- 
cesses in the Quantum Chromodynamics, and in a more or less detailed presentation can 
be found in any textbook. Less widely known, is the fact that from the theoretical point 
of view this picture is not complete. An indication that some contributions may be miss- 
ing, comes from the asymptotic nature of the perturbative series in (0). This series is 
non-Borel-summable, which means that any attempt to attribute a quantitative mean- 
ing to the sum of the series in (|3]) would produce an exponentially small imaginary part 
~ zexp{— const • tt /a s (Q 2 )}, which is to be cancelled by the imaginary part coming from 
nonperturbative contributions. Thus, separation between perturbative and nonperturba- 
tive pieces in the cross section as the ones which contribute to the coefficient function 
and the parton distribution, respectively, cannot be rigorous. A modern discussion of the 
asymptotical properties of the perturbation series in QCD can be found in [[[], 0. 

In addition to imaginary exponential corrections which must cancel identically against 
the corresponding nonperturbative contributions, the coefficient functions may acquire 
also real exponential corrections, which potentially produce observable effects. In this 
talk I shall report on recent results by I. Balitsky and myself || indicating that these 
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corrections are indeed present. We have found that the deep inelastic cross section indeed 
possesses exponential contributions of the form F(x) exp[— 4irS(x)/a s (Q 2 )}, where S(x), 
F(x) are certain functions of Bjorken x, which we are able to calculate in a certain 
kinematical domain. Since the experimental data are becoming more and more precise, it 
is of acute interest to find a boundary for a possible accuracy of the perturbative approach, 
which is set by nonperturbative effects. Our study has been fuelled by recent findings 
of an enhancement of instanton-induced effects at high energies in a related problem of 
the violation of baryon number in the electroweak theory 0, §J. In the case of QCD the 
instanton-induced effects could turn out to be numerically large at high energies, despite 
the fact that they correspond formally to contributions of a very high fractional twist 
exp(-47rS(x)/« s (g 2 )) ~ (A 2 QCD /Q 2 ) bS ^ . 

To be precise, in this talk I consider the contribution to the structure functions com- 
ing from the instanton-antiinstanton pair. Contributions of single instantons are only 
present for higher-twist terms in the light-cone expansion and are less interesting. A 
Ringwald-type enhancement || of the instanton-induced cross sections at high energy 
can compensate the extra semiclassical suppression factor exp(— 2-k/q) accompanying 
instanton-antiinstanton contributions compared to single-instanton ones. In such case 
the II terms become the leading ones owing to a bigger power of the coupling in the 
preexponent. 

As it is well known, the instanton contributions in QCD are in general infrared- 
unstable. In a typical situation integrations over the instanton size are strongly IR- 
divergent. Our starting point is the observation that this problem does not affect cal- 
culation of instanton contributions to the coefficient functions. Let us introduce for a 
moment an explicit IR cutoff Am to regularize the integrals over the instanton size. Then 
the contribution of the instanton-antiinstanton pair to the cross section can be written 
schematically as 

Q 2 a(Q 2 ) ~ (Am/A QCD ) 2b + {A QCD /Q) 2bS{x) . (4) 

The second term in (|]) gives an IR-protected contribution. It depends in a nontrivial 
way on the external large momentum and is identified unambiguously with a contribution 
to the coefficient function. The first term contributes to the parton distribution. To be 
precise, one should separate in the first term the contributions coming from instanton 
sizes above and below the reference scale /i, and to add the contribution of small-size 
instantons to the coefficient function. Schematically, one has in this way 



(Air/Aqcd) = (/V Aqcd) + \ {Am/ Aqcd) - ^/A QCD r\ . (5) 

However, this reshuffling of the Q 2 -independent contribution between the coefficient and 
the parton distribution does not affect the observable cross section. It is analogous to 
an ambiguity in the separation between contributions to the coefficient function and to 
the parton distribution in perturbation theory, induced by possibility to use different 
regularization schemes (e.g. MS instead of MS, etc.). Hence, we can concentrate on 
contributions of the second type in (f|), which are IR-protected. 

2. The distinction between the instanton-induced contributions to the coefficient 
functions, which are given by convergent integrals over the instanton size, and the con- 
tributions to parton distributions, given in general by IR-divergent integrals, becomes 
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Figure 1: The contribution of the instanton-antiinstanton pair to the cross section of hard 
gluon-gluon scattering (a), structure function of a gluon (b) and of a quark (c,d). Wavy lines 
are (nonperturbative) gluons. Solid lines are quark zero modes in the case that they are ending 
at the instanton (antiinstanton), and quark propagators at the // background otherwise. 

especially transparent for an example of the cross section of hard gluon scattering from a 
real gluon, see Fig. la, considered in detail in |3]]. This process is not directly physically 
relevant (e.g., the cross section is not gauge invariant), but it serves as a good toy model. 
Following Zakharov J7[, we calculate the cross section by means of the optical theorem. 
The trick is to evaluate the contribution to the functional integral coming from the vicin- 
ity of the instanton-antiinstanton configuration in Euclidian space, and calculate the cross 
section by the analytical continuation to Minkowski space and by taking the imaginary 
part. Each hard gluon is substituted by the Fourier transform of the instanton field in 
the singular gauge at large momentum, and brings in the factor || 



All dependence on the hard scale comes in this case through the explicit power of Q 2 in 
front of the divergent integral. This is a typical contribution to the parton distribution 
- in the present case to the probability to find a hard gluon within a soft gluon. The 
second term gives rise to a completely different behavior. The cross section is given in this 
case by the following integral over the common scale of the instanton and antiinstanton 
Pi ~ pi and over their separation R in the cm. frame || |3|: 





The first term in (^) produces a power-like divergent integral over the instanton size p. 




(7) 
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Three important ingredients in this expression are: the factor exp(— 2Qp), which comes 
from the two hard gluon fields, the factor exp(ERo), which is obtained from the standard 
exponential factor exp[i(p + q)R],E 2 — (p + q) 2 by the rotation to Minkowski space, 
cf. 0, and the action S{^) evaluated on the instanton-antiinstanton configuration. The 
normalization is such that S"(£) = 1 for an infinitely separated instanton and antiinstanton, 
and £ is the conformal parameter [TIJ 



P1P2 

Writing the action as a function of £ ensures that the interaction between instantons is 
small in two different limits: for a widely separated // pair, and for a small instanton put 
inside a big (anti)instanton, which are related to each other by the conformal transforma- 
tion. In the limit of large £ the expansion of >S(£) for the dominating maximum attractive 
// orientation reads JT0[ 




S{t)= l-^ + O(m(0/£ 4 ) (10) 



where the l/£ 2 term corresponds to a slightly corrected dipole-dipole interaction. Thus, 
the action S(£) decreases with the distance between instantons, so that the instanton and 
the instanton effectively attract each other. This attraction results in the exponential 
increase of the cross section — the effect found by Ringwald [H]. Further terms in the 



expansion of the action can be obtained by the so-called valley method |TT| , and a typical 
solution (conformal valley) |10| gives a monotonous function of the conformal parameter, 
which turns to zero at R — > 0. In the traditional language, the valley approach corre- 
sponds to the summation of all so called "soft-soft" corrections arising from the particle 
interaction in the final state. Main problem is in the evaluation of "hard-hard" corrections 



121 , which come from particle interaction in the initial state. These corrections are likely 
to decrease the cross section, and in physical terms must take into account an (exponen- 
tially small) overlap between the initial state, which involves a few hard quanta, with the 



semiclassical final state ||13|| . Thus, the instanton-antiinstanton action is substituted by 
an effective "holy grail" function, which determines the leading exponential factor for the 
semiclassical production at high energies, and which received a lot of attention in recent 



years. Unitarity arguments [|14l O suggest that the decrease of the action will stop at 



values of order S(£) ~ 0.5. In a recent preprint [16[ Diakonov and Petrov argue that 
S(£) indeed decreases up to the value 1/2 at a certain energy of order of the sphaleron 
mass, and then starts to increase, so that the semiclassical production cross section is 
resonance-like. The question seems to us to be not settled finally. In this study, we have 
taken the value S = 1/2 as a reasonable guess for the residual suppression, and assumed 
that the behavior of the "true" function S(£) for S(£) > 1/2 is close to that given by the 
conformal valley [PJ. The latter assumption is supported by numerical studies, e.g. in 
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To the semiclassical accuracy the integral in (||) is evaluated by a saddle-point method. 
The saddle-point equations take the form || 

Aix 

Qp* = — r ^(£.-2)S'&) + &S , (&). 
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Figure 2: The non-perturbative scale in deep inelastic scattering (instanton size 1 ), 
corresponding to the solution of saddle-point equations in ([ll]) as a function of Q and for 
-0.5-0.6 (£*~3-4). 

8tt 



where S"(0 is the derivative of £(£) over £, and p*,£* are the saddle-point values for the 
instanton size and the conformal parameter, respectively. 

Neglecting in (O) the terms proportional to b = (ll/3)iV c — (2/3)n/, which come from 
the differentiation of the running coupling (and produce a small correction), one finds 

R 2 l + x 

pi i-x 

Qp* = -^-^ (12) 

A numerical solution of the saddle-point equations in (|TlD for the particular expression of 
the action S(£) corresponding to the conformal instanton-antiinstanton valley is shown 
in Fig. 2. Note that the difference between the hard scale Q 2 and the effective scale for 
nonperturbative effects p~ 2 is numerically very large. This is a new situation compared 
to calculations of instanton-induced contributions to two-point correlation functions, see 



e.g. |g, [17], |18|| , where the size of the instanton is of order of the large virtuality. The effect 
is that the instanton-induced contributions to deep inelastic scattering may turn out to 
be non-negligible at the values Q 2 ~ lOOOGeV 2 , which are conventionally considered as a 
safe domain for perturbative QCD. 

In the case of hard gluon-gluon scattering it is easy to collect all the preexponential 
factors (to the semiclassical accuracy). The result for the scattering of a transversely 
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polarized hard gluon from a soft gluon reads 



2E 2 a± 
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-d 2 



[1 - x) 2 + x 2 
x 2 (l — x) 



-It 



13/2 



' 27T 



21/2 



exp 



47T 



+ 2b S&) 



(13) 



It is expressed in terms of the saddle-point values of p and £. Here d 0.00363 (for 
rif = 3) is a constant which enters the expression for the instanton density 



d 



-Cxesxp[nfC 3 -N c C 2 )], 



(14) 



d = 0.466, C 2 = 1.54, C 3 = 0.153 in the MS scheme. 

Note that the preexponential factor in ( [L3|) is calculated to the leading semiclassical 
accuracy, i.e. in the limit x — > 1, while in the exponential factor exp{ — (47i/a s + 2b)S(^)} 
we used the expression for S(£) suggested by the valley method, which is likely to be 
a good approximation provided S(£) > 1/2, alias x > 0.25 — 0.3, cf. (]T2]). With this 
restriction, and at a s (p*) ~ 0.3 — 0.4, the expression on the r.h.s. of (|13D reaches values 
of order 10~ 2 - 10°, which means that at Q 2 ~ 100 - 1000GW 2 and x < 0.25 - 0.40 the 
nonperturbative contribution appears to be significant. 

3. Similar contributions are present in the structure functions of deep inelastic 
lepton-hadron scattering, but the calculation turns out to be much more involved 0. 
The situation proves to be somewhat simpler for the case of deep inelastic scattering from 
a real gluon. To this purpose we need to evaluate 



T, 



i / d x e^(A%p),X\T{j^x)U0)}\A a (p),\) 



W, 



— ImT, 



7T 



(15) 



q 2 



F L (x,Q 2 ) + 



PuPv ViAv + %Pv 



pq 



pq 



Q 



+ 9^^) F 2 (x,Q' 



Q 



There are two technical problems to be solved. First of all, the separation of the finite 
contribution to the coefficient function under the background of a divergent contribution 
to the parton distribution is no longer given by a simple formula in (^]). Instead, we 
extract the contribution of interest by making an analytic continuation of integrals over 
the instanton size p from negative values of p 2 . For a typical integral we write down, e.g. 



dp [ 



T(A) 



(T 2 + p 2 ) 



,2\A 



2i sin[7r(A — p — 

p 2 + ie 
T 2 + p 2 + ie / 



n)} 

A 



— c.c. 



dp {-p 
T(A 



2\li+n— A— 1 



p — n)T(p + n) 



(T 



(16) 



The second, and main complication, comes from the necessity to consider the quark 
propagator in the II background, see the diagram in Fig. lb. Neglecting corrections of 
order p 2 /R 2 in the preexponential factor, we can make use of the cluster expansion ||, 
and keep the first nontrivial term only: 

d_ 



{x\V7?V lT \0) 



dz (^V^Vil^^-^lVsV^lO). 
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Here and below the subscript '1' refers to the antiinstanton with the size p\ and the 
position of the center xj = R + T, and the subscript '2' refers to the instanton with 
the size p2 and the center at xj = T. We use conventional notations V = V M cr M and 
V = V^cr^, etc, where a° a = (— ia, 1), a^a — (+^0", 1), an d 0" are the standard Pauli 
matrices. The expressions for quark propagators at the one-instanton (antiinstanton) 
background are given in [|TJJ . 



We have found that to the leading accuracy in the strong coupling the instanton 

contribution in Fig. lb comes from the following integration regions in the coordinate 
space: 

2 / \2 2 

z ~ (z — X) ~ X 

(x-R-T) 2 + p\ ~ T 2 + p 2 2 ~x 2 

(z — R — T) 2 + p\ ~ (z - Tf + p 2 ~ x 2 /a s 

T 2 ~ # 2 ~ _ p 2 _ ^2^ ^ 17 ^ 

We remind that all the calculation is done in Euclidian space, and the evaluation of 
integrals by means of the analytical continuation effectively corresponds to the integration 
over negative values of p 2 , see flT6|) . 

Hence the integration over z in the cluster expansion of the quark propagator can be 
done in the "light-cone" approximation: 

F(z) ir 2 /-i , F(jx) _ 



(a;-z) 4 z 4 x 4 Jo 7(1-7) 

where F(z) is an arbitrary function containing all other possible denominators like (z — 
R — T) 2 + p\ etc. This is a major simplification compared to the general case. 

After a considerable algebra, we obtain the following answer for the II contribution 
to the structure function of a real gluon: 

Fi '{x,Q ) = 



q 9x 2 bSfo)[bS(&)-l] 

( 4tt 
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(19) 



where the expressions for p* and ^* coincide to the ones given in (O) . To our accuracy, we 
find that the instanton- induced contributions obey the Callan-Gross relation ^(x, Q 2 ) = 
2xF!(s,Q 2 ). 

The expression in (|1^) presents our main result. It gives the exponential correction to 
the coefficient function in front of the gluon distribution of the leading twist in @. The 
exponential factor is exact to the accuracy of (0). The preexponential factor is calculated 
to leading accuracy in the strong coupling and up to corrections of order 0(1 — x). The 
corresponding contribution to the structure function of the nucleon is obtained in a usual 
way, making a convolution of (|T9"D with a distribution of gluons in the proton at the scale 

Pi 
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In the language of the operator product expansion, the answer in flT3 ) can be rewritten 
as the exponential nonperturbative contribution to the high n ~ Tr/a s moment of the 
structure function 

Mf\n,Q 2 )= C dxx n Ff\x,Q 2 ) (20) 
Jo 

Taking n = h- 47r / 'o. s {p*) with n of order unity, one can evaluate the integral in ( p0| ) with 

F (G) 

(x, Q 2 ) given in (|19D by the saddle point method. The saddle-point equation reads, 
approximately 

S%)(^ 2 -4)/4 (21) 



n 



where is related to the saddle-point value of x by (|T2"D. The lowest possible value ~ 3, 
corresponding within the valley-method approximation to S(£*) ~ 1/2, 5"(£*) ~ 0.27, thus 
yields n ~ 0.3. Hence the expression in (|Tj^) is applicable to the evaluation of instanton- 
induced corrections to the coefficient functions in front of local operators starting from 
n ~ 10- 12. 

Leading contribution to ( |T9| ) in the perturbation theory is due to the mixing with the 
flavor-singlet quark distribution and is given by the box graph: 



Ipert 



E2r 2 , -2i a s(Q , , 
e q [x + x I — z In 



2vr 



Q 2 x 

fJ, 2 X 



(22) 



where in order to compare to the instanton contribution in ( |i~9|) one should choose the 
scale p to be of order p*. 

The instanton-antiinstanton contribution to the structure function of a quark contains 
a similar contribution shown in Fig.lc. The answer reads 



Fi q \x,Q 2 ) 



E4 + 



128 



d 2 7T 9 / 2 



16 



81xHS(Q[bS(&-l] \e 



x 



2tt 



15/2 



exp 



47T 



+ 26 5(6 



(23) 



However, in this case additional contributions exist of the type shown in Fig. Id. They 
are finite (the integral over instanton size is cut off at p 2 ~ x 2 /a s ), but the relevant 
instanton-antiinstanton separation R is small, of order p. This probably means that the 
structure of nonperturbative contributions to quark distributions is more complicated. 
This question is under study. The answer given in (|23|) presents the contribution of the 
particular saddle point in (|T2|) . 



4. The instanton-induced contribution to the structure function of a gluon in ( |I~9D 
is shown as a function of Bjorken x for different values of Q ~ 10 — lOOGeV in Fig. 3. The 



contribution of the box graph in (22) is shown by dots for comparison. The low boundary 
for possible values of Q is determined by the condition that the effective instanton size 
is not too large. At Q — lOGeV we find p* ~ lGeV^ 1 , cf. Fig. 2. This value is 
sufficiently small, so that instantons are not distorted too strongly by large-scale vacuum 
fluctuations. Another limitation is that the valley approach to the calculation of the "holy 
grail" function S(£,) is likely to be justified at S'(^) > 1/2, which translates to the condition 
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Figure 3: Nonperturbative contribution to the structure function F\(x, Q 2 ) of a real gluon ( |l9l ) 
as a function of x for different values of Q (solid curves). The leading perturbative contribution 



(22) is shown for comparison by dots. The dashed curves show lines with the constant effective 



value of the action on the II configuration. 

that x > 0.3 — 0.35. Numerical results are strongly sensitive to the particular value of the 
QCD scale parameter. We use the two-loop expression for the coupling with three active 
flavors, and the value = 290MeV which corresponds to = 240MeV Q. The 
corresponding value of the coupling at the scale of r-lepton mass is a s {m T ) = 0.29. The 
new ALEPH data [F2TJ indicate a somewhat higher value a s (m T ) = 0.33 ± 0.05. Since the 



dependence on the coupling is exponential, the 20% increase of a s (p*) induces the increase 
of the cross section by an order of magnitude! Together with uncertainties in the function 
S(£) and in the preexponential factor, this indicates that the particular curves given in 
Fig. 3 should not be taken too seriously, and rather give a target for further theoretical 
(and experimental?) studies to shoot at. 

To summarize, we have found that instantons produce a well-defined and calculable 
contribution to the cross section of deep inelastic scattering for sufficiently large values of 
x and large Q 2 ~ 100 — lOOOGeV 2 , which turns out, however, to be rather small — of order 
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1CT 2 — 1CT 5 compared to the perturbative cross section. This means that the accuracy 
of standard perturbative analysis is sufficiently high, and that there is not much hope 
to observe the instanton-induced contributions to the total deep inelastic cross section 
experimentally. However, instantons are likely to produce events with a very specific 
structure of the final state, and such peculiarities may be subject to experimental search. 
The dominating Feynman diagrams in our calculation correspond to 27r/a(p*) ~ 15 gluons 
and 2nf — 1 = 5 quarks in the final state, with the low virtuality of order p~ x ~ 1 GeV, 
and which are produced in the S-wave in the cm. frame of the colliding partons. It 
is not likely that these can be resolved as minijets, and we rather expect a spherically 
symmetric production of final state hadrons. The effect is likely to be resonance-like, that 
is present in a narrow interval of values of Bjorken x of order 0.25-0.35 (in the parton- 
parton collision). In any case, finding of an instanton-induced particle production at high 
energies is a challenging problem, and further theoretical efforts are needed to put it as a 
practical proposal to experimentalists. 

5. I would like to thank Ian Balitsky for a very rewarding collaboration. It is a 
pleasure to acknowledge also useful discussions with M. Beneke, M.A. Shifman and V.I. 
Zakharov. 
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